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An exact spin mapping is identified to simplify the recently proposed hard-core boson description 
(Demler and Das Sarma, Phys. Rev. Lett., to be published) of the bilayer quantum Hall system at 
filling factor 2. The effective spin model describes an easy-plane ferromagnet subject to an external 
Zeeman field. The phase diagram of this effective model is determined exactly and found to agree 
with the approximate calculation of Demler and Das Sarma, while the Goldstone-mode spectrum, 
order parameter stiffness and Kosterlitz-Thouless temperature in the canted antiferromagnetic phase 
are computed approximately. 
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Quantum Hall systems have attracted prolonged in- 
terest, because they not only exhibit the extremely in- 
teresting quantized Hall effect, but also prove to be the 
ideal arenas to study many other important subjects, 
such as Anderson localization, Luttinger liquid behav- 
ior, itinerant magnetism, to name a few. Recently, much 
attention has been, focusing on multicomponent quan- 
tum Hall systems!!! Such components may be the spins 
of electrons which are not frozen out when the external 
magnetic field is not too strong, or layer indices in multi- 
layered system. Novel physics such as even-denominator 
fractional quantum Hall statesjjtl spin-ferromagnetism 
and topological excitations ,□ and spontaneous inter-layer 
coherencaJH have been found. 

Very recently, a novel canted antiferromagnetic phase 
has been proposed theoretically in bilayer quantum Hall 
systems at filling factor M—,= 2, based primarily on 
Hartree-Fock calculations. Elo There is some experimen- 
tal support for the existence of such a phase.tLil Later, 
an interesting attemptLHl was made to go beyond the 
Hartree-Fock theory. In this attempt, one first truncate 
the Hilbert space in such a way that the sum of the num- 
ber of electrons for a given Landau orbital is always two 
for both layers, so that the total charge degree of freedom 
is not allowed to fluctuate. The problem is further simpli- 
fied by keeping only two lowest energy states for a given 
orbital. With these simplifications and mapping the re- 
duced Hamiltonian onto hard-core bosons on a lattice, 
Demler and Das Sarmali3 were able to obtain a phase dia- 
gram that agrees well with Hartree-Fock in the absence of 
disorder, and predict some novel properties of the system 
with disorder. The advantage of this formalism is that 
certain local correlations are treated non-perturbatively, 
and effects of disorder are easy to incorporate. Thus it is 
complementary to Hartree-Fock. 

The purpose of the present paper is to show that the 
hard-core bosons introduced in Ref. [ [l2| form an exact 
Schwinger boson representation of an effective spin-1/2 
system. In terms of the spin-1/2 operators, the reduced 
Hamiltonian of Ref. [ [l2| takes a particularly simple 
form; it describes an easy-plane ferromagnet subject to a 



magnetic field along the z direction, while the XY plane 
is the easy plane. The advantage of this spin-mapping 
is that a lot of insight from, and theoretical techniques 
developed for spin models (which are not available in the 
original hard-core representation) can be applied to the 
present problem. In the spin model, it is clear that the 
phase diagram obtained in Ref. [ ^] based on approx- 
imate variational wave functions that only satisfy the 
hard-core condition in average, is actually exact, within 
the reduced Hamiltonian. More importantly, the stan- 
dard spin- wave theory applied to the present model yields 
the exact collective mode spectrum for the singlet and 
triplet phases (again within the reduced Hamiltonian), 
as well as the approximate Goldstone mode spectrum 
and order parameter stiffness in the canted phase. From 
the order parameter stiffness we can also estimate the 
Kosterlitz-Thouless temperature for the canted phase. 

We begin with the reduced hard-core boson Hamilto- 
nian, Eq. (4) of Ref. [pf: 



- -(cosfl + sm9) 2 J2(4^h + t)v 3 v\u). 
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Here i is the index of a complete and orthogonal set of 
localized intra-Landau level orbitals which form a square 
lattice. t\ and v\ are boson creation operators which 
create a pair of electrons in orbital i; the pair created by 
t\ have their spins polarized along the external magnetic 
field to take advantage the Zeeman energy, while the pair 
created by v\ occupy the symmetric combination of the 
two states in two layers to minimize the tunneling energy, 
and form a spin singlet. These are two low-energy local 
states being kept. Since the number of electrons is fixed 
to be 2 for each orbital, they must satisfy the hard-core 
condition 
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where is the Zeeman splitting, e is the dielectric 
constant of the system, I is the magnetic length, and 
e c = {e 2 /e£)y^/2(l - e~* x l £ ' ' Erf[d/(V2£)]) is the local 
charging energy, in which d is the layer spacing and Erf 
is the error function. 

The constraint Eq. (|2|) implies that the t and v hard 
core bosons form a Schwinger boson representation^ of 
a spin-1/2 operator, for each lattice site. We may define 
the spin-up state to be the state in which a t boson is 
present, and the spin-down state to be the one with v 
boson present. The spin-1/2 operators are related to the 
hard-core boson operators in the following way: 
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It is important to remember that these effective spin- 
1 /2 degrees of freedom correspond to the singlet and spin- 
polarized electron pair states, and are not directly related 
to the physical spins of the individual electrons. In terms 
of these effective spin-1/2 operators, the reduced Hamil- 
tonian Eq. dll) takes a particularly simple form: 
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in which 
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where Z is the coordination number of the lattice and 
for a square lattice, Z = 4. Since Jii > J± > 0, Eq. 
( |l0| ) describes an easy-plane (XY plane) ferromagnet, 
in the presence of a Zeeman field along the z direction, 
B z . Again, B z is an effective field that couples to these 
effective spins, which should not be confused with the 
physical magnetic field. 

The existence of three phases of the original system 
is particularly easy to understand in terms of the spin 
model ([H]). When B z is large and positive, all the ef- 
fective spins are pointing up in the ground state (|vl/|) = 
| || • • ■ I)), which means all the pairs are in the t state. 
This corresponds to the fully spin-polarized ferromag- 
netic (FPF) phase. When B z is large and negative, all 
the effective spins are pointing down (|*x) = | H ••• |)), 



which means all the pairs are in the v state, this corre- 
sponds to the spin singlet (SS) phase. In between, how- 
ever, the effective spins would like to develop a sponta- 
neous polarization in the XY plane, to take advantage of 
the J|i coupling (since it is stronger than J±). This state 
with a spontaneously broken XY or U(l) symmetry cor- 
responds to the canted antiferromagnetic (CAF) phase. 
These three phases are schematically illustrated in Fig. 
1. 

The phase diagram of ( |l0| ) may be determined exactly, 
because the spin polarized (in either up or down direc- 
tion) states are exact eigen states of the Hamiltonian, 
due to the fact that Sf ot = J2i Sf is a conserved quan- 
tity. So are the single magnon excitations on top of them. 
When l^i) is the ground state, the exact single magnon 
spectrum is 

E(k) — B z + 2J± — Jii (cos k x a + cos k y a), (14) 



where a-.= yl-nt is the lattice constant for the square 
lattice.O This FPF phase becomes unstable against the 
CAF phase when the spectrum becomes soft at k = 0, 
namely B z + 2J± - 2J« = 0, or 



E v -E t + J(l - sin26») = 0. 



(15) 



This determines the exact phase boundary separating the 
FPF and CAF phases, within the model (JTTj). The phase 
boundary separating the CAF and SS phases may be 
obtained similarly: 



E v -E t + J(sin 9 + cos 9) 2 = 0. 



(16) 
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FIG. 1. Configurations of effective spins in the three 
phases. All spins are pointing up in the fully spin-polarized 
(FPF), and pointing down in the spin singlet (SS) phase. In 
the canted antiferromagnetic (CAF) phase, the spin polariza- 
tion is tilted away from the z direction by an angle <f>, and a 
component in the XY plane is spontaneously developed. 
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These agree with the results of Ref. [ [L2|, although 
in Ref. [ |l2|] it was not clear the phase boundaries ob- 
tained were exact within the reduced Hamiltonian, as 
variational wave functions that only satisfy the hard-core 
condition on average were used. 

In the CAF phase, it is no longer possible to obtain the 
exact ground state or magnon excitation spectrum, due 
to the presence of quantum fluctuations. In this case, we 
may perform an approximate spin- wave analysis. As a 
first step, consider the following variational wave func- 
tion, that describes the state in which all the effective 
spins are still fully polarized, with the polarization tilted 
away from the z direction by an angle <f> (see Fig. 1): 
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This variational wave function (with <j> being the varia- 
tional parameter) is different from that used in Ref. [l2| 
and has the advantage that the hard-core condition is 
satisfied exactly, in the original hard-core boson repre- 
sentation. The optimal 4>, which minimizes ((f>\H\ip), is 
found to satisfy 



sin(b(B z - 2AJ cos cb) = 0, 



(18) 



in which A J = Jy - J± = (J/2) sin(20). The FPF and 
SS phases correspond to the trivial solutions <f> = and 
(j) = 7T, while in the CAF phase we have 



= cos- 1 [B z /(2AJ)}. 



(19) 



The phase boundaries can also be determined from the 
relative stability of these solutions, which agree with the 
above spectrum analysis. 

In order to perform a spin- wave analysis, we need to 
rotate the spin quantization axis so that it is parallel to 
the spin polarization of \<p), after which the Hamiltonian 
reads 

H = -B z ^(cos 4>St + sin <j>Sf) ~ J\\ ^ S, • S^- 

m 

+ AJ^cos^S? + emjSf) ■ (cos<pS* + sin<£SJ). (20) 

(ij) 

Now use r- the standard Holstein-Primakoff 
transformationEj to map spin operators to boson oper- 
ators: 

S$ = S 71 i, 

Sr = V2Sf s ( ni )al 

S+ = sf2Sf s { ni )ai, (21) 

where a\ is the boson creation operator for site i, rij = 
a\a.i, S = 1/2 is the size of the effective spin, and 



fs(n) 
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(22) 



Keep terms that are linear or quadratic in boson opera- 
tors in H, we obtain the linearized spin wave Hamiltonian 
(up to a constant): 



Hlsw = 2J|| ^2 ojoj + - AJsin 2 0^(a|a] + a^a,) 
« to) 

- - ^AJsin 2 0) X>J aj + a)oi). (23) 

(ij) 

The terms that are linear in a cancel due to the condi- 
tion Eq. ( |l9| ) . This boson Hamiltonian may be diagonal- 
ized by first going to momentum space and then perform 
a Bogliubov transformation. The diagonalizcd Hamilto- 
nian takes (up to a constant) the following form: 



Hlsw — ^ E^b^b^, 



(24) 



where the Bogliubov boson b is related to a by 

Ok = coshak&k + sinhak&Lk! (25) 

in which 



— y 



(26) 



and the angle «k satisfies 

tanh2a k = -A k /e k , (27) 

where 

Ak = [(A J/2) sin 2 <f>](cos k x a + cos k y a), 
ek = 2J|| — [J|| — (A J/2) sin 2 <^>](cos k x a + cos k y a). (28) 

The spin wave spectrum is 

£k=V^-A£. (29) 



Taking the long wave length limit k — > 0, we find a linear 
spectrum: — c|k| expected for a Goldstone mode, 
with the Goldstone mode (spin wave) velocity 

J a sin ( 



c = a sin 4>\ I Ju A J = 



A/sin 26>(1 + sin 26»). (30) 



Quantum fluctuations of spin-waves reduces the mag- 
netization in the CAF phase. Within the linear spin-wave 
theory, the magnetization per unit cell at T = is 

m =- — — V(4«k) = i - — y 1 



k Vi - (A k A k ) 

(31) 

The susceptibility along z direction is 
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(32) 



in which we used Eq. (|19|). 

The important quantity, order parameter stiffness p s , 
is reiated to Xz and c viaE£l 



2 / 2 

X^c /a 



(33) 
and 



It is straightforward to calculate p s numerically, 
some numerical data are presented in Fig. 2. 

Assuming the short-distance physics at finite temper- 
ature to be similar_to that of the classical XY model 
on a square lattice,cl the Kosterlitz-Thouless transition 
temperature may be estimated to be 



T, 



KT 



0.9 Ps 



(34) 



Thus in these systems, the typical Tkt is of order 
TfCT ~ 0.2J - 0.005f^ - 0.5K. This is qualitatively 
consistent, although somewhat smaller, than the Hartree- 
Fock estimated The fact that the present estimate gives 
a smaller number may be due to fluctuations effects that 
have been taken into account in the reduced Hamiltonian 
(which were left out in Hartree-Fock). 
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FIG. 2. p 3 (in unit of J) as a function of (j>, for several 
different 9's. From bottom up, 9 = n/32,n/16, 7r/8, 7r/6, 7r/4. 
p s is invariant under 9 =>• n/2 — 6, or <fi 7T — <j>. 

We emphasize that charge fluctuations have been left 
out in the reduced effective spin Hamiltonian studied 
here, and it describes neutral excitations of the system 
only. In particular, the skyrmions and vortices of the ef- 
fective spin model are charge-neutral. This is very differ- 
ent from the bilayer systems at filling factor 1, where all 
vortices must carry a half-integer chargeBEZI On the other 
hand, the present system can support charged skyrmions 
and vortices, which are being studied .ETlj 
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